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Abstract 
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Nomenclatiirp 

nondlmensional elastic axis 
location, positive aft of 
midchord 

airfoil semi -chord 
nondimensional aileron hinge 
line location 

nondimensional damping matrix 
control distribution vector 
plunging d.o.f., positive 
downward from elastic axis 
displacement, velocity, and 
acceleration control gains 
respectively ’ 

nondimensional stiffness matrix 
^ss of the airfoil per unit sjan 

free stream Mach number 
nondimensional mass matrix 
nondimensional sensor location, 
positive aft of midchord 
aerodynamic load vector 
radii of gyration of 
airfoil about clastic axis and of 
aileron about hinge axis, 
respectively 

o+iw, Laplace transform variable 
time; w^jt, nondimensional time 
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tree stream velocity; U/bov,, 
nondimensional flight speed 
nondimensional distances 
from elastic axis to airfoil mass 
center and hinge line to aileron 
mass center, respectively 

displacement and state vectors, 
respectively cci-ors, 

pitching d.o.f., positive 
leading edge up pw>M.ive 

tJ]?nnn'’lS‘''’^23 positive 

trailing edge down 

control surface conmand 

control surface viscous 

damping ratio 

m/»pb2. airfoil mass ratio 
n/b. nondimensional plunging 

eia°;ti’c aMs^^'^ 
nondimensional sensor 
plunging displacement, 
focatiw **ownward from sensor 

free stream air density 
uncoupled natural 

abo3t^Al^!^°^ plunging, pitching 

DitcMnl ^ aileron 

pitching abut hinge axis. 

respectively 

Introduction 
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present stupy. furtPer results’ aT p'reLtL “r 
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the MBBA-3 airfoil. Three degrees of freedom 
10 . 0 . f.) are considered; plunge, oitch anH 
pitch. Response analyses are performed 

eauaU^]‘^or“lV- ^"^®9'‘2t1n9 the structural 
equations of motion along with the unsteadv 

of transonic code LTRAN2- 
similar Identification technique 

amHoH 1.°^ Bennett and DesmaraisH ^is 

identify* thp® time-marching response curves to 
loentify the aeroelastic modes. Estimated 

inT^V^r? "“«<es are plotted 

in a root -locus type format in the comolex 

a statp’c analyses are performed using 

a state-space aeroelastic model, termed the Fade 

approximants of the 
unstea^ aerodynamic forces. These forces are 
^f^'-’bed by an interpolating function 

Thp^p^na^” harmonic transonic aerodynamic data. 
The Padd model is written as a set of linear 

eouation^®''* coefficient, differential’ 

equations. These equations are solved in the 

the yielding eigenvalues compared in 

estimates’ Vhp aching modal 

estimates. The objective of the stability 

Inri'ifu^i'* investigate the applicability of 
locally linear aeroelastic modeling to airfoils 

stahVm®°'”'^ '^eP'^esentative Padd model 

stability results for the NACA 64A010 airfoil 

Padd model may 

alternatively be solved in the time domain with 
appropriate initial conditions yielding the 
aeroelastic displacement time-histories. 13^ 

Open-loop stability and response anaivses 
are performed to determine the behavior of'^the 

^ function of flight 
speed. Time-marching response calculations are 
performed at three different flight speeds to 
investipte subcritical, critical, and super- 
critical flutter conditions. Frequency and 
aeroelastic modes identified from 
2del Sits" compared with Padd 

;>,-o 1-1°sed-loop stability and response analyses 
plwtfP to investigate application ofS 

analvses^To time-marching 

rt suppression using active 

is intentionally 

effects P“'‘P°ses. Aeroelastic 

effects due to simple, constant gain, feedback 
control laws utilizing displacement, velocity 

variety of control gains. ^ 

. . The objective of this study is to gain 

insight into the control of airf^oil stabi1i-ty 
and response in transonic flow. The focus of 
the present effort is on time-marching transient 
response solutions with aeroelastic Mdal 
identification. The calculations presented here 
reveal some interesting aeroelastic behavior 

svsterwHh The aeroelastic 

elaSc Itahn'^V''® controls. Transonic aero- 
nf V^n response characteristics 

are rade^”'*°’ ^ discussed and comparisons 


Time-Marching Response Analysis 

Response analyses are performed by 
simultaneously integrating the structural equa- 


tions of motion along with the unstp.rfu 

transonic code 

[M](X} + [C]{xj + [K]jx} = {p} + {Gjgj. ( 1 ) 
to nondimensional time u^t. respect 

•*c = Xd«s * Ms ^ Ms (2) 

vSlocit?’ and displacement. 

resSSCelv^ £ acceleration control gains, 
Inn - The sensor measured oluna- 

thl airVon «« Plaled along 

rfoil chord to obtain a measure of airfoil 
plunge and pitch motions from 


«s = iMjlX} 


(3) 


jrs isif ' 

EQ"^ 2 ^ ° control law. 

2si cho’rd ^ SdPi'ace of 

aerodvnan?iV f ®°"Trol surface 

aerodynamic forces are utilized to alleviate 
flutter instability. Details of the contmi 

system and equations may be found in Ref. 12. 

expressing the control law Eq. (2) in 

terms of the airfoil motion, the aeroelastic 

motion may be written in general 
form for time-i ntegrat ion as 

[M*]lX) + [C*](X) + [K*]jx} = jp) ( 4 , 

where [M*] = [m] - K^{g}iHj ( 5 a, 

(c*] = [c] - KylGjiHj (55, 

[■^*] = [Kj - Kp{G}iHj ( 5 J., 

Details of the time-marching response solution 
procedures were given in Ref. 13. 

A modal identification technique is used to 

oSseTf Jho frequency, amplitude, aSd 

phase of the aeroelastic modes from the time- 

displacement response histories. A 

De la?ai Bennett anJ 

fit th. ; ' 'east-sci.ares c.r.e- 

E;«to's !" hrS“ 

m (— ) .T 


^(T) = a + j; e “o J [a. cos (— , . f 
j=l J “a'j 

* '’j T J 

rt 


(6) 
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The damping and frequency of the complex modes 
thus obtained 


(2_). + i . (i_\. 

“a J “a J “a J 


(7) 


are estimates of the aeroelastic eigenvalues and 
mcy then be directly compared in the complex 
s-plane with those computed by the Padd model. 

As an example. Fig. 1 shows a typical 
response analysis of the MBB A-3 airfoil at 
M = 0.765 and Cj^ = 0.58. The LTRAN2-NLR time- 
marching pitching response a and the modal curve 
fit using Eq. (6) are shown in the top part of 
the figure. Only the data used for curve 
fitting (144 points) is shown. In the lower 
part of the figure are the two component aero- 
elastic modes, identified from the modal fit. 


Pad6 Model Stability Analysis 

Stability analyses are performed using a 
linear eigenvalue analysis of the Padd model. 
This model was formulated by curve-fitting the 
unsteady aerodynamic forces by a Padd approxi- 
mating function. These approximating 
functions are then expressed as linear differen- 
tial equations which, when coupled to the struc- 
tural equations of motion, lead to the first- 
order matrix equation 


U1 = I- |Z} = [A]{Z) + {B}b ( 8) 

where {z} is the state vector containing 
displacement, velocity, and augmented states; 
s/o)a is the complex eigenvalue. By expressing 
the control law Eq. (2) in terms of the state- 
vector (Z), Eq. (8) is easily solved using 
standard eigenvalue solution techniques. 
Details of the Padd model formulation are given 
in Ref. 12. 


Results and Discussion 


Transonic aeroelastic stability and 
response analyses were reported In Ref. 9 for 
three airfoil configurations: NACA 64A006, NACA 

64A010, and MBB A-3. All three airfoils are 
among those proposed by AGARD for aeroelastic 
applications of transonic unsteady aerody- 
namics. Aeroelastic behavior of the MBB A-3 
airfoil was studied at the design Mach number 
M = 0.765 and at (1) zero mean angle of attack 
a = (P; and (2) the design steady lift coeffi- 
cient Cl = 0.58. The mean angle of attack 
TOcessary to match Cj^ = 0.58 using LTRAN2-NLR 
is a = 0.860. In this report, representative 
time-marching response histories and s-plane 
stability root-loci for the MBB A-3 airfoil at 
M = 0.765 and c^ = 0.58 are presented. Aero- 
elastic parameter values selected are the same 
as those used in Refs. 12 and 13 and are listed 
in Table 1. 


Table 1. Aeroelastic parameter values for 
stability and response analyses. 


(ji./iii ~ 0.3 

n a 

X - 0.2 
a 


= 0.5 

M = 50.0 

Xp = 0.008 

ah = “0.2 

r^ = 0.06 

C3 " 



Steady pressure solutions are required as 
aerodynamic initial conditions for unsteady cal- 
culations. As an example, the MBB A-3 steady 
pressure distribution computed using LTRAN2-NLR 
at M 0,765 and c^ = 0-58 is shown in Fig. 2 
along with the airfoil contour. This steady 
pressure data compares well with the experimen- 
tal results of Bucciantini et al.l6 There is 
a relatively weak shock wave on the airfoil 
upper surface near 55% to 60% chord. Steady 
pressure distributions for all of the cases 
considered indicate shock locations in the range 
of 50% to 65% chord and shock strengths that are 
well within the range of applicability of 
transonic small -disturbance theory. 9 


Open-Loop 

Open-loop stability and response analyses 
are performed to determine the behavior of the 
aeroelastic modes as a function of flight speed. 
These calculations were first performed using 
the Padd aeroelastic model with results plotted 
in a flight speed root-locus format. As an 
example, open-loop root-loci as a function of 
flight speed are presented in Fig. 3 for the MBB 
A-3 airfoil. With increasing flight speed, the 
torsion branch moves to the left in the stable 
left-half of the complex s-plane. The aileron 
branch is also stable throughout, moving up and 
to the left as flight speed is increased. The 
bending dominated root-locus becomes the flutter 
mode at a nondimensional flutter speed of 
Up = 2.370. In addition, the Padd model also 
predicts the divergence speed, where an aerody- 
namic lag root moves onto the positive real 
axis, as Up = 3.320. This divergence pheno- 
menon Is similar to that reported by 
Edwards, where static divergence of a 
typical section in incompressible flow occurred 
due to the emergence of a real positive root 
from the complex s-plane origin. In Ref. 14, 
this root appeared in addition to the original 
structural poles and was also predicted using a 
Padd model. Table 2 compares the flutter and 

divergence speeds with those computed using 
linear aerodynamic theory. Flutter speeds from 
a p-k method flutter analysis are also tabulated 
for further comparison. These p-k method 
calculations were performed to assess the 

accuracy of the Padd model flutter solution. 
Padd model flutter speeds compare well with the 
p-k method values. 
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Table 2. Suirmary of open-loop flutter and 

divergence speeds for the MBB A-3 
airfoil at M = 0.765. 


AIRFOIL 

FLUTTER 

METHOD 

FLUTTER 

SPEED 

U* 

DIVERGENCE 

SPEED 

^D 

"d 

MBB A-3 
sO.58 

Fade 

2.370 

3.320 

1.401 


P-k 

2.403 



LINEAR 

THEORY 

Padd 

2.711 

3.611 

1.332 

p-k 

2.729 




Time-marching transient responses were then 
obtal ned at three d1 f f erent f 1 1 ght speeds to 
investigate subcritical, critical, and super- 
critical flutter conditions. These calculations 
were performed independent of the Padd approxi- 
mation, primarily for verification purposes. 
Transient response histories are shown In Fig. 4 
for three values of flight speed near the 
flutter speed. A one-percent semi-chord plunge 
displacement was used as initial condition. 
Time-responses for U/Up = 1.0 In Fig. 4 are 
nearly neutrally stable corresponding to the 
flutter condition. The flutter speed used for 
time-marching analyses was the p-k method 
value. Aeroelastic transients for 

U/Up = 0.844 and U/Up = 1.156 in Fig. 4 show 
converging (subcritical ) and diverging (super- 
critical) oscillatory behavior, respectively. 

Damping and frequency estimates of the 
aeroelastic modes were determined from the time- 
response histories of Fig. 4. These estimates 
are compared with Pad^ model eigenvalues for the 
bending and torsion modes in Fig. 5. In 
general, the Pad6 results compare well with the 
time- marching modal frequency and damping 

values. Agreement for the lower frequency bend- 
ing mode is, in general, better than that for 
the higher frequency torsion mode. Frequency 

and damping estimates for the aileron mode are 
much less accurate than those for bending and 
torsion. The above discrepancies are attributed 
to: (1) the Pad^ model is less accurate at the 
higher damping ratios involved; (2) the 
structural integration is performed with fewer 
time-steps per cycle of the higher frequency 
modes; and (3) since the aileron mode has much 
more damping in comparison with bending and 
torsion, its contribution to the responses dies 
out quickly leaving very little information for 
modal curve-fitting. 


Closed-Loop 

The Pad6 aeroelastic model was used to 
study the effects of active feedback control on 
the aeroelastic modes at the flight speed equal 
to the open-loop flutter value. These calcula- 
tions were performed primarily to determine the 
range and sign of the control gains required to 
stabilize the flutter mode. Padd stability cal- 
culations were also performed at the same Mach 
numbers using linear subsonic aerodynamic theory 
for comparison with transonic Padd model 
results. Effects due to displacement, velocity, 


and acceleration feedback were further assessed 
using time-marching response analyses for 
comparison with Padd results. Finally, an 
example is presented illustrating the 
effectiveness of velocity feedback upon flutter 
suppression at speeds above and below the open- 
loop flutter speed. 


Effects of Active Control on Flutter Mode .- 
Closed-loop analyses are first performed at the 
f 1 1 ght speed equal to the open-1 oop f 1 utter 
value, to determine the behavior of the aero- 
elastic modes as a function of control gain. 
Effects of active control on the bending 
dominated flutter mode are of primary interest. 
Control gain root-loci computed using the Padd 
model for displacement, velocity, and accelera- 
tion control laws are shown in the left, center, 
and right of Fig. 6, respectively. Correspond- 
ing root-loci computed using linear subsonic 
theory are given in Fig. 7. In each case there 
are three open-loop aeroelastic roots or "poles" 
corresponding to the three structural d.o.f. *s. 
The two poles in the stable left-half of the 
complex s-plane correspond to the torsion and 
aileron modes; the pole that lies on the 
imaginary axis is the bending dominated flutter 
mode. These aeroelastic roots move toward 
"zeros" of the aeroelastic transfer function as 
the control gains Kp, K\j, or are 

monotonically increased or decreased. Solid- 
line root-loci indicate the sign of the control 
gai n that stabilizes the flutter mode whereas 
the dashed-line loci indicate the opposite 
effect. Comparison of LTRAN2-NLR results of 
Fig. 6 with the linear theory results of Fig. 7 
indicates that frequency and damping values are 
significantly different for transonic and 
subsonic cases, although the overall root-locus 
trends are similar. This lack of agreement 
between linear theory results and LTRAN2-NLR 
root-loci illustrates the importance of 
including transonic effects in aeroelastic 
stability calculations. 

Displacement feedback with negative control 
gains easily stabilized the flutter mode as 
shown in Figs. 6 and 7. Negative displacement 
feedback, however, results in static divergence 
of the same nature found in the open-loop cases. 
For the MBB A-3 airfoil at M = 0.765 and 
Cl = 0.58, this occurs for Kq £ -0.5. 

Divergence root locations are shown in Figs. 6 
and 7 for K[) = -1.0. Also, positi ve 

displacement feedback destabilized the flutter 
mode. 

Velocity feedback with positive control 
gains stabilized the bending dominated flutter 
pole and increased the damping of the torsion 
mode. The transonic results of Fig. 6 indicate 
that there is more damping in bending and 
aileron modes, and much less in torsion than 
predicted by linear subsonic theory in Fig. 7. 
Also, negative velocity feedback destabilized 
the flutter mode. 

Acceleration feedback with positive control 
gains stabilized the flutter pole and decreased 
damping and frequency of the torsion mode. 
Also, negative acceleration feedback destabil- 
ized the flutter mode. 
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Effects Due to Displacement Feedback . - 
Aeroelastic effects due to displacement feed- 
back, ^s» were further studied by obtaining 
time-marching transient responses for 
Kq = -1*0 and detailed Pad6 stability results 
for the range of displacement gains 
-1.0 <. Kq <. These calculations were 

performed prTmarily to investigate the diver- 
gence phenomenon predicted by a positive real 
aerodynamic lag root in the Padd model. The 
three d.o.f. time-histories computed using 
LTRAN2-NLR at U/Up = 1.0 and Kq = -1.0 are 
given in Fig. 8. These responses show that 
flutter has been suppressed and that the aero- 
elastic displacements are Indeed divergent. 
Comparisons of time-marching frequency and 
damping estimates with Padd eigenvalues are 
shown in Fig. 9. Good agreement is found in 

predicting this aeroelastic divergence. In 
addition, negative displacement feedback 
increased the damping in bending and decreased 
the torsion mode frequency. The lower frequency 
bending mode damping and frequency values for 
Kd = -1.0 compare well; results from the time- 
marching modal fits for the higher frequency 
torsion and aileron modes were not reliably 
determined. 

Padd stability calculations were then 
performed for the range of gains 
-1.0 ^ Kq < 0.0 to reveal the effects of nega- 
tive dispracement feedback on flutter and 
divergence speeds. These more detailed results 
are shown in Fig. 10 for the MBB A-3 airfoil. 
Negative displacement feedback slowly increases 
the flutter speed and first rapidly increases 
the divergence speed. For gains Kq between 
approximately -0.2 and -0.5, divergence is 
eliminated. In this region, the aerodynamic lag 
root previously causing divergence no longer 
crosses through the origin. For Kn ^ -0.5, 
divergence reappears as the critical mode of 
instability at speeds much lower than flutter. 
The flutter speed is increased approximately 14% 
for Kp =* -0.5. 


Effects Due to Velocity Feedback . 
Aeroelastic effects due to velocity feedback , 
^s» were further investigated by obtaining 
transient responses for successively increased 
velocity feedback. These calculations are first 
performed at the flight speed U* = Uj?. Values 
selected for the control gains were Ky = 3.0, 

6.0, and 9.0, the same gains shown by the tri- 
angles in the center portion of Fig. 6. 
LTRAN2-NLR response results are shown in Fig. 
11. These response histories show that as Ky 
is increased the dominant motion is consistently 
more damped and of higher frequency. Also, 
positive velocity feedback significantly 
increased the nonrat ional part^^ of the plunge 
displacement responses. Here, the ^-transients 
oscillate about an asymptotically decaying 
rather than a zero mean value. 

Damping and frequency of the aeroelastic 
modes determined from the LTRAN2-NLR displace- 
ment transients are plotted in Fig. 12 along 
with Pad^ model results. The flutter mode 
damping and frequency are both successively 
increased for values of Ky = 3.0, 6.0, and 

9.0. The torsion mode damping is significantly 


Increased with little effect upon the torsion 
mode frequency. Comparison of Padd eigenvalues 
with time-marching modal estimates shows general 
agreement. Discrepancies between the two sets 
of results are attributed to the same reasons as 
in the open-loop case. Also, the torsion mode 
for Kv = 6.0 and 9.0 could not be identified 
from the time-marching transients because of the 
large increase in damping. 

Velocity feedback with positive control 
gains suppressed flutter for the range of flight 
speeds investigated. Some sensitivity with 
respect to mean angle of attack has been 
observed though in calculations using velocity 
feedback. 12 por the MBB A-3 airfoil at 
M = 0.765 and zero mean angle of attack, the 
flutter speed was increased by only 11% for the 
gain Ky = 9.0. Additional calculations also 
showed that velocity feedback does not affect 
the divergence speed. Therefore for the 
Cl = 0.58 case, the aeroelastic instability 
becomes that of static divergence. Since the 
control law with velocity feedback was found to 
be generally the most effective in raising 
flutter speeds, i2 more detailed calculations 
were performed. These results are described in 
a following section entitled “Illustrative 
Example. " 


Effects Due to Acceleration Feedback . - 
Aeroelastic effects due to federation feed- 
back, were further investigated by obtain- 
ing transient responses for successively 
increased acceleration feedback. These calcula- 
tions are performed at the flight speed 
U* = Up. Values selected for the control 
gains were K/\ = 6.0, 12.0, and 18.0, the same 
gains shown by the triangles in the right por- 
tion of Fig. 6. LTRAN2-NLR response results are 
shown in Fig. 13. These response histories show 
that as K/\ is increased the dominant motion is 
consistently more damped and of lower frequency. 
Also, a higher frequency transient becomes more 
visible in the a and e responses for the larger 
values of This is due to the decreased 
damping of the torsion mode as shown in the 
frequency and damping comparisons of Fig. 14. 
Positive acceleration feedback increased damping 
in bending, decreased damping in torsion, and 
lowered the frequencies of both modes. Compari- 
son of Pad6 eigenvalues with time-marching modal 
estimates shows general agreement. Results for 
the lower frequency bending mode are in better 
agreement than results for the higher frequency 
torsion mode. 

Acceleration feedback with positive control 
gains raised the flutter speed approximately 25% 
for the maximum gain studied, K^ = 18.0. 
Additional calculations showed that acceleration 
feedback does not affect the divergence speed. 


Illustrative Example . - Additional aero- 
elastic analyses were performed at speeds above 
and below the open-loop flutter speed using the 
velocity feedback control law since this control 
law was found to be generally the most 
effect! ve. These results are presented to 
further illustrate application of the Padd 
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and time-niarching methods. Padd flight speed 
root- loci are shown In Fig. 15 for values of the 
velocity gain Ky = 3.0, 6.0, and 9.0. Posi- 
tive velocity feedback significantly increased 
the damping of the torsion root-locus and 
decreased the damping of the aileron mode. The 
bending root-locus for Ky = 3.0 reflects an 
increase in flutter speed of approximately 29% 
(Up = 3.050), while for Ky = 6.0 and 9.0 
flutter has been eliminated. Since velocity 
feedback does not affect the divergence speed, 
the aeroelastic instability becomes that of 
static divergence at the uncontrolled divergence 
speed = 3.320. This divergence speed is 
approximately 40% greater than the open-loop 
flutter speed as shown in Table 2. 

Time marching response calculations were 
then performed for U/Up = 1.266 using LTRAN2- 
NLR to complement and verify the Padd model 
results. These response histories shown in 
Fig. 16, correspond to the U* = 3.0 diamond- 
symbol eigensolutions of Fig. 15. These aero- 
elastic transients are bending dominated and 
converging for all three values of the control 
gain Ky. As Ky is increased, the responses 
become more stable indicating a monotonic 
Increase in bending mode damping, at approxi- 
mately the same frequency. Again, positive 
velocity feedback significantly increased the 
nonrational part of the plunge displacement 
responses. 

Effects of velocity feedback on bending and 
torsion modes at U/Up = 1.266 are shown in 
Fig. 17. Because of the large damping increase 
in torsion due to Ky and the relatively high 
frequency of the aileron mode, only the bending 
mode was determined from the aeroelastic tran- 
sients of Fig. 16. Comparison of the two sets 
of bending root-locus results in Fig. 17 shows 
very good agreement. Positive velocity feedback 
increased bending mode damping at approximately 
the same frequency. 

Finally, open and closed-loop flight speed 
root-loci are shown in Fig. 18 for both bending 
and torsion modes. The velocity feedback gain 
used to generate the closed-loop results is 
Ky = 9.0, the same value used in the right 
portion of Fig. 15. Pade model and time- 
marching solutions for the bending mode compare 
well. The control gain Ky = 9.0 increased the 
torsion mode damping for the range of flight 
speeds 0.844 ^ U/Up jC 1.266, such that the 
results do not permit comparison between Padd 
and time-marching solutions. Fig. 18 

illustrates again that positive velocity 
feedback has eliminated the flutter instability 
for the MBB A-3 airfoil at M = 0.765 and 
Cjt = 0.58. 


Conclusions 

Transonic aeroelastic stability and 
response calculations have been performed for 
2-D airfoils using the LTRAN2-NLR transonic 
small -disturbance code. Representative results 
are presented here for the MBB A-3 airfoil at 
the design condition M = 0.765 and Cj^ = 0.58. 
Emphasis is placed on time-marching transient 
response solutions with modal identification for 


comparison with Padd model eigenvalues. Exten- 
sion of the time-marching transient response 
analysis to a three degree of freedom aero- 
elastic system including active controls is 
demonstrated. Accurate frequency and damping 
estimates were obtained using a complex 
exponential least squares curve-fit of the aero- 
elastic responses. Good agreement was found 
between s-plane eigenvalues calculated using the 
Padd aeroelastic model and time-marching modal 
estimates. Therefore, locally linear aero- 
elastic modeling was found to be applicable to 
2-D airfoils in small -disturbance transonic 
flow, for cases that are within the range of 
validity of the transonic codes. The importance 
of including transonic effects in aeroelastic 
stability calculations was illustrated by the 
differences between root-loci calculated using 
linear subsonic theory and LTRAN2-NLR. Although 
the overall root-locus trends are similar, 
frequency and damping values are significantly 
different for transonic and subsonic cases. 

Effects due to simple, constant gain, 
control laws utilizing displacement, velocity, 
or acceleration sensing were studied. Displace- 
ment feedback with negative control gains 
suppressed flutter but caused airfoil responses 
to diverge. Comparison of Pade eigenvalues with 
time-marching frequency and damping estimates 
showed good agreement in predicting this aero- 
elastic divergence. In addition, negative dis- 
placement feedback increased the damping in 
bending, and decreased the torsion mode 
frequency. Detailed Padd stability calcula- 
tions revealed the effects of negative displace- 
ment feedback on flutter and divergence speeds. 
Velocity feedback with positive control gains 
suppressed flutter. Padd stability results 
compared well with time-marching frequency and 
damping values. In general, positive velocity 
feedback Increased damping in both bending and 
torsion modes, and significantly increased the 
nonrational part of the plunge displacement 
responses. Here, the plunge transients 
oscillated about an asymptotically decaying 
rather than a zero mean value. Acceleration 
feedback with positive control gains also 
increased flutter speeds. Again, Padd results 
compared well with time-marching modal 
estimates. In general, positive acceleration 
feedback increased damping in bending, decreased 
damping in torsion, and lowered the frequencies 
of both modes. 
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Fig. 1 Modal curve fit of aeroelastic 

time-marching pitching response of the 
MBB A-3 airfoil at the open-loop flutter 
speed U/Up = 1.0, K/v = 18.0, 

H = 0.765, and C£ = 0.58. 



Fig. 2 Steady pressure distributions for the 
MBB A-3 airfoil at M = 0.765 and 
Cjt = 0.58. 
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Fig. 3 MBB A-3 open-loop flight speed root-loci 
at M » 0.765 and Cj^ = 0.58. 



Fig. 4 Effect of flight speed on MBB A-3 

open-loop time-marching displacement 
responses at M = 0.765 and * 0.58. 



Fig. 5 MBB A-3 flight speed root-loci for 

bending and torsion modes at M = 0.765 
and Cl “ 0.58. 



Fig. 6 MBB A-3 control gain root-loci at the 
open-loop flutter speed (U* = Up), 

M = 0.765 and Cj^ = 0.58. 



Fig. 7 Control gain root-loci at the open-loop 
flutter speed = Uf) computed using 
linear subsonic aerodynamic theory at 
M = 0.765. 
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/3*io axio 



F1g. 10 Effect of negative displacement 
, .. X ji. 1 , UDD feedback on MSB A-3 flutter and 

F1g« 8 Effect of displacerient feedback on MBB divergence speeds at M ■ 0»765 and 

A-3 time-marching displacement responses . , 5,58, 

at the open-loop flutter speed * 

U/Up “ 1.0, M ■ 0.765 and ■ 0.58. 




% 



Fig. 9 Effect of displacement feedback on MBB 
A-3 bending and torsion inodes at the 
open-loop flutter speed U/Up ** 1.0, 

M « 0.765, and c^ * 0.58. 



TIME (W^t) IN 27T RADIANS 


Fig. 11 Effect of velocity feedback on MBB A-3 
time-marching displacement responses at 
the open-loop flutter speed U/Up = 1.0, 
M » 0.765 and Cj^ » 0.58. 
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Fig. 12 Effect of velocity feedback on MBB A-3 
bending and torsion inodes at the 
open-loop flutter speed U/Up ■ 1.0, 

M » 0.765 and Cji « 0.58. 


Fig. 14 Effect of acceleration feedback on MBB 
A-3 bending and torsion modes at the 
open-loop flutter speed U/Up - 1.0, 

M * 0.765 and Ci « 0.58. 




Fig. 15 MBB A-3 velocity feedback, flight speed 
root-loci at M ® 0.765 and cj^ = 0.58. 


Fig. 13 Effect of acceleration feedback on MBB 
A-3 time-marching displacement 
responses at the open-loop flutter 
speed U/Up « 1.0, M « 0.765 and 
Cjf, * 0.58. 
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Fig. 18 MBB A-3 open and closed (Ky » 9.0) loop 
flight speed root-loci for bending and 
torsion modes at M * 0.765 and Cj^ “ 0.58. 


Fig. 16 Effect of velocity feedback on MBB A-3 
time-marching displacement responses at 
U/Up « 1.266, M « 0.765 and c^ « 0.58. 





Fig. 17 Effect of velocity feedback on MBB A-3 
bending and torsion modes at 
U/Up » 1.266, M « 0.765 and « 0.58. 
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